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Abstract. 

We give a formulation of quantum cosmology with a pressureless dust 
and arbitrary additional matter fields. The system has the property that its 
Hamiltonian constraint is linear in the dust momentum. This feature provides a 
natural time gauge, leading to a physical hamiltonian that is not a square root. 
Quantization leads to Schrodinger equation for which unitary evolution is directly 
linked to geodesic completeness. Our approach simplifies the analysis of both 
Wheeler-deWitt and loop quantum cosmology (LQC) models, and significantly 
broadens the applicability of the latter. This is demonstrated for arbitrary scalar 
field potential and cosmological constant in LQC. 



PACS numbers: 98.80Qc, 04.60Kz, 04.60Pp 

1. Introduction 

In attempts to formulate a non-perturbative theory of quantum gravity one of 
the obstacles is the so-called problem of time. Its manifestation in the canonical 
formulation of general relativity is that the Hamiltonian vanishes, which is a direct 
consequence of the time reparamctrization invariance of the theory. For matter with 
the usual kinetic energy, the Hamiltonian constraint is quadratic in all momenta, 
including the pure gravity term. 

There appears to be an obvious way to address this problem: impose an explicit 
time gauge fixing to obtain a true Hamiltonian that describes the evolution of physical 
degrees of freedom. This explicit deparametrization brings with it some difficulties: 
the true Hamiltonian is generically a square root, which is difficult to define as an 
operator. A second problem is that different time gauge fixings give radically different 
true Hamiltonians, so in the final analysis the quantum theory ends up being time- 
gauge dependent. 

In a recent paper [T] we presented a complete non-perturbative quantization of 
gravity that solves both problems. This approach uses a single timelike dust field 
(motivated by the Brown and Kuchaf (BK) form) 2 coupled to general relativity and 
standard model matter. Like the BK action, the canonical theory has the feature that 
the Hamiltonian constraint is linear in the dust momentum. However, unlike the BK 
action, it also has the remarkable feature that in a natural time gauge fixing the true 
Hamiltonian is not a square root. Furthermore, this Hamiltonian is identical in form 
to what would be the Hamiltonian constraint without the dust field. It is this feature, 
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combined with the well-developed kinematical framework for loop quantum gravity, 
that leads to a complete quantum theory. 

In this paper we apply this development to quantum cosmology in both the 
Wheeler-DeWitt and LQG quantization. We show that (i) the true Hamiltonian is 
first order, and in the LQC case it poses no self-adjoint extension problem, (ii) the 
non-zero cosmological constant is analytically solvable for both signs, and (iii) any 
type of matter is included. 

Our approach differs significantly from existing LQC models, all of which use the 
scalar field as a clock; this relies on the scalar field momentum being a constant of 
motion, and leads to a Hamiltonian constraint operator that is second order in time 
(and so is not a time-dependent Schrodinger equation). For these reasons it is difficult 
to implement scalar field time for arbitrary scalar potential, including polynomial or 
slow roll cases used in inflationary models. 



2. General relativity with dust 

We consider the class of the flat isotropic spacetimes whose metric is of the form 

ds 2 = -N 2 dt 2 + a 2 (i)(dx 2 + dy 2 + dz 2 ) . (1) 

The theory is given by the Einstein-Hilbert action with timelike dust field T and an 
arbitrary matter lagrangian C m (assumed to be at most second order in time) : 

1 

4G 



S = I d A x\/ —gR — I d 4 .T\/— g C v 



1 / - (2) 

- 2 / d X V=gM(g ab d a Td b T + 1) . 

where M is a Lagrange multiplier. The last term resembles the Brown-Kuchaf action 
[2], but contains only one scalar. 

The canonical formulation of this action gives a Hamiltonian constraint linear in 
the dust momentum px conjugate to T, together with the usual spatial diffeomorphism 
constraint [T| . The linearity in px makes the theory resemble a parametrized system, 
and suggests the gauge T = t. In this gauge the true Hamiltonian is —pr, given by 

H = - PT = H G + H m (3) 

where Hq and H m are the gravitational and matter parts of the Hamiltonian 
constraint coming from the first two terms in the action. In this natural gauge, 
quantization gives the time dependent Schrodinger equation 

= [H G + H m )y (4) 

Since Hq and H rn are manifestly independent of the dust variables, the physical 
hamiltonian is time independent. The spatial diffeomorphism constraint of course 
vanishes identically in the homogeneous sector. 

The dust effectively provides a cosmic time, therefore there is a direct relation 
between unitarity and geodesic completeness of the resulting spacetime in this time 
parameter; if classical singularities are not removed by quantum dynamics, loss of 
unitarity occurs precisely at the singular points. In the following we will see illustrative 
examples of this. 
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3. Quantum cosmology with dust 

As a first example we consider FRW cosmology with only dust, ie. H m = in (j4|), and 
give a quantization of this system using both the Wheeler-DeWitt (WDW) and the 
background independent LQG method. These methods differ in the choice of Hilbert 
space used. This will demonstrate a first advantage over scalar field time used in all 
previous LQC work: classically the singularity is reached in finite dust time, so it 
allows a direct comparison of the mechanism of singularity resolution in background 
independent quantization, and the lack of it in WDW quantization. 

The phase space variables are the canonical pair {a,P a } — 1 where a is the scale 
factor Jl}. Classically a is positive, therefore to have phase space topology IR. 2 ^] one 
needs to equip a with an orientation. To do this we can either work with the original 
variables or use a canonical transformation to the oriented volume v (— a~ 1 a 3 ) l§| and 
its conjugate momentum 6, which satisfy the Poisson bracket {v, b} = 2. In these 
variables H G = (3irG/2a)b 2 \v\. 



3.1. Wheeler-DeWitt quantization 

In WDW quantization one uses the standard Schrodinger representation. The 
canonical variables (v, b) are promoted to operators defined on the Schwartz space 
<S(R) € Hg = L 2 (M., dv), where s denotes restriction to the functions symmetric in v. 
Since the classical Hamiltonian is linear in w, a naive realization of the Hamiltonian 
operator in momentum space is 

H G := Sini^oT 1 b 2 d b . (5) 

It acts on a dense domain Hg of functions spanned by the eigenfunctions 
9(±x) exp(ikx), where k G R and x = —1/6. This operator however breaks the 
positive definiteness of the classical Hamiltonian since its spectrum Sp(Hc) = R- An 
obvious way to fix this problem is just to restrict attention to the positive part of ((5|). 
However there is a rigorous way to arrive at this conclusion in the following steps: (i) 
Begin with the the natural factor ordering [3j 0] 

H G = ^VW\b 2 VW\- (6) 

This is a symmetric positive definite operator on <S(M), (ii) analysis of its deficiency 
subspaces shows that it admits a family of extensions labeled by e I/3 € U(l), and 
that the spectrum of each extension Hp with domain T-Lp is nondegenerate. (hi) Due 
to this fact, for each extension Hp, there is an invertible map P_p : Hp — > Hp to a 
certain auxiliary space such that Hg is mapped onto the positive part [if g] + of Hg- 
(The auxiliary spaces Hp are determined by the selfadjoint extensions Hp of Hg, 
and construction of the map uses the fact that the symmetric state is encoded in one 
orientation 6(v)^f(v).) 

In each extension the physical states are 

\PpV](x) = J dk^(k) [6(x)e lkx + 6{-x)e l P e lkx ] , (7) 

| This gives a simple transformation between the configuration and momentum representations. 
§ a s; 1.35£pj, This choice synchronizes the variable v with what is used in LQC. a is expressed in 
terms of the Barbero-Immirzi parameter and the LQC area gap [3]. 
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where $ G L 2 (R, fcdfc), with the inner product 

(Pp*\Pp$) = I dx\P^¥](x)[id x ][P^](x) . (8) 
Jr 

Time evolution is given by 

tf t (fc) = e l "( fc )(*-*°)^ o (fc) , u (k) = 37r4iO -1 * , (9) 

This is a plane wave packet moving with the constant speed v := Sirlp^a^ 1 till it hits 
x = 0, where it gets rotated by an extension dependent phase e t/3 . 

The dynamics of the system is described by the observable V — \a 3 \. Under the 
mapping V — P^VPp 1 , it takes a simple form involving only d x and x 2 I. Its time 
evolution and dispersion is analytically determined: [jj] 

(V){t) = V(t) = 2a- 1 (k) [Z-Kll x {t - t s )] 2 + 2aa 2 x , k = kl , (10) 

where t s corresponds to the (always existing) point where (:id x x:) = (the colons 
denote symmetrization) and a 2 can be related with the dispersion of x at this moment. 
In the same way one can determine the time dependence of Hubble parameter H . 
Physically these trajectories follow the classical ones, and describe a contracting 
universe until the big crunch singularity at the time t s . At that moment V becomes 
comparable to its dispersion. In this sense the singularity V — is reached dynamically. 

To determine the evolution past the singular point one needs to know the 
selfadjoint extension. This is equivalent to imposing additional boundary conditions 
at the singularity x = 0; it cannot be determined by an analysis of the state for t < t s . 
Given a choice of extension it is possible to continue the evolution. The resulting 
trajectory for t > t s describes a universe that expands out of the singularity and 
subsequently follows a classical trajectory. Thus we see that in WDW quantization 
the big bang singularity is dynamically reached, and evolution beyond it is not 
uniquely determined without additional boundary data. In this sense the the WDW 
quantization does not resolve the singularity. 



3.2. Loop Quantum Cosmology 

We now consider the same problem in LQC. To quantize the system we follow the 
improved dynamics prescription of [3 ]. As in the WDW case, we can apply directly 
the results of kinematical quantization [3l [5]. The canonical variables are again 
(v, b), but the physical Hilbert space is now H = L 2 (R, d/is), where K is the Bohr 
compactification of the reals and d/x# is the Haar measure on it. The (physical) 
gravitational Hamiltonian Hq is a second order difference operator in v. With the 
same factor ordering as in (|6]), it takes the form 

Hg = -^VWlW-^'fVWl , M\v) = \v + 1). (11) 
8a 

This operator is non-negative definite and, unlike in WDW quantization, it is 
essentially selfadjoint [6]. The fact, that Hg does not couple the sets v < 0, v = 
and v > permits with construction of the mapping P : He — > Ha into an auxiliary 
space analogous to that for WDW quantization. Due to the selfadjointness of Hq its 
construction here is even simpler. The auxiliary Hamiltonian now takes the form 

H G := PHgP- 1 = [Zin^a" 1 sm 2 (b)d b } + , (12) 



| The explicit dependence of the first term on a is compensated by analogous dependence in (k). 
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where due to discreteness in v the momentum b is periodic with period n. The auxiliary 
Hilbert space T-Lq consists of states 

[P#](x) = f dk^(k)e- tkx , f eL 2 (I+,Hi), (13) 

where now x = — cot (6). The inner product (in new coordinate) is given by the 
same formula ©. At this point it is worth noting that Hq is a 1st order differential 
operator so its spectrum is nondegenerate. This fact is crucial in identifying the 
images of the mapping of the original Hilbert space basis to the auxiliary one. In 
the models with a scalar field an analogous construction would lead to the 2nd order 
operator with degenerate spectrum. This make the identification of the bases much 
more involved and required an alternative approach [5] (which cannot be extended to 
more complicated systems.) 

The physical evolution is again given by but due to selfadjointness the 
evolution across x = is now unique. In particular the expectation value of the 
volume is now 

V(t) = 2a- 1 (k) [iTi(l x {t-t s )] 2 + 2ad 2 x + 2a{k) , (14) 

where the meaning of t s and a\ is the same as in (fTUl) . We note here the presence of 
a nontrivial minimum volume V m = 2a(k), which implies dynamical resolution of the 
singularity through the bounce, as in the models with scalar field time [31 [S]. 

In summary we have seen that the dust reference frame gives exactly solvable 
models in WDW and LQC, a feature present in earlier models [5]. We now turn to 
more complicated models, where the real advantage of dust time becomes apparent. 



3.3. LQC with a cosmological constant 

FRW cosmologies with non- vanishing cosmological constant and scalar field time have 
been investigated in detail in LQC [71 18]. In these works the appearance of the 
square root in the evolution operator rendered the models analytically unsolvable. 
Furthermore, for the case A > the Hamiltonian was not essentially selfadjoint; the 
multitude of extensions complicated the analysis even further, in particular raising an 
issue of dependence of global evolution on the choice of time [4]. As we will see, the 
situation simplifies considerably with dust time. 

With the cosmological constant the gravitational Hamiltonian is 

H A = H G + ^A\v\ , (15) 
16a 

(where p c rj 0.81ppi is the critical energy density of LQC [3].) The polymer 
quantization of this operator is straightforward [7j : It is positive definite for A < 
and essentially selfadjoint for all values of A. Remarkably, the cosmological constant 
term does not disrupt the properties required for the mapping Pa : Hg — > Ha between 
the physical and auxiliary Hilbert spaces. The auxiliary Hamiltonian is 



H A := PaHaPK 1 = H g - ^Ad b , (16) 



This is simple enough to determine its eigenfunctions analytically. For all values of A 
these can be expressed in the form e lkx ( b ) ^ where x(b) depends analytically on A. 
A < 0: For this case we have 



x(b) = arctan 



^/(3/9 c |A| + 247r4i)/(3/> c |A|) tan(6) 



(17) 
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This inherits the periodicity properties of b. Therefore the spectrum of the auxiliary 
Hamiltonian [-Ha] + is discrete. Nondegeneracy of its spectrum allows again to 
relate the original Hilbert space basis elements with the auxiliary counterparts in 
a straightforward way. The (auxiliary) physical states are 

Ha 3 [Pa*] (*) = ]T * n e 2mx , (18) 

and the time evolution is given by 

*«(*) = *„(Oe^*-*°) , cu n = £ v /3 /9c |A|(3 Pc |A| + 24^2 i) . (19) 

The standard observables used for A = also map in a straightforward way into an 
auxiliary system. In particular V :— PaVP^ 1 of V acts directly on the spectral profile 
in terms of the shift and multiplication operators. 

The calculation of expectation values reproduces the dynamics qualitatively 
similar to the one of the scalar field system: the universe goes through an infinite 
chain of quantum bounces and classical recollapses. Now however, unlike in [7] due to 
exact uniformity of Sp(i? a) the evolution is exactly periodic, thus the wave packets 
do not spread out between the cycles of the evolution. 

A > 0: As in the model with scalar field we observe qualitative changes at the critical 
value A c = 8TrGp c . The most interesting case is for < A < A c . Let us write 
A as A = A c sin 2 (/3), with (3 G [0,7r/2]. Then the domain of b ([0,7r]) admits two 
uncoupled deSitter sectors. Furthermore the spectrum of auxiliary Hamiltonian (|TB)) 
is degenerate, which makes the identification of the basis vectors of Ha considerably 
more involved. Fortunately a careful analysis of the asymptotic properties of the 
eigenfunctions of H\ allows them to be uniquely identified. This gives that Sp(if a) = 
M with the eigenfunctions 

e +|fc ,(6) = 0(sin(6)-sin(/3)) e -+( b ) , 

e_ ]fc |(6) = 0(sin(/3) - sm(b))e ix -^ , 
where the coordinates x± are related to b by 

tan(/3)[tanh(x±)] Tl = tan(6) , x± G R . (21) 

We have thus two sectors, one of positive and the other of negative energy, supported 
respectively on sin(6) > sin(/3) and sin(fe) < sin(/3). In each sector, physical 
states are of a form analogous to (fT5|) . with k G ifl, a basis defined via (|2U)) 
and * G L 2 (M., \k\dk) . They evolve as free plane wave packets in the respective 
coordinates x±. Time evolution is again given by eqn. Q , now with the frequencies 
uj(k) — Siri^a^ 1 cos 2 (/3)fc. The volume operator is in this case involves d x± and 
multiplications by hyperbolic functions of x± . One can thus apply the methods of [5] 
to calculate its expectation value and dispersion analytically. 

The positive energy sector gives analytical results similar to those found with 
scalar field time in LQC, ie. a contracting deSitter-like universe bouncing once into an 
expanding one. In addition we obtain negative energy solutions which give DeSitter 
universes with "phantom" dust and classical bounce, i.e. due to classical negative 
energy. The last feature has interesting implications for the case A > A c : there is a 
large solution space of phantom dust solutions with classical bounce (A = A c ) or cyclic 
evolutions consisting of an infinite chain of classical bounces and quantum recollapses 
(A > A c ). These cases however are just mathematical curiosity and are not of physical 
interest. 

% Each sector is represented by k > and k < respectively. 
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4. Other applications 

As we have seen, the use of dust time considerably simplifies the analysis of systems 
with cosmological constant. Its usefulness is not however restricted to only these: the 
physical Hamiltonian is self-adjoint for any non- exotic matter, including Yang-Mills 
gauge fields, and scalar fields with any potential term. This follows from the fact that 
Hg is self-adjoint for any value of the cosmological constant j9], and that the only 
gravitational variables the matter terms contain are factors of the volume operator 
and its inverse. This holds also for other spatial topologies. 

Infiaton potential: One of potentially important applications is the analysis of the 
cosmologies with massive scalar field, including those with infiaton potential. In this 
situation the phase space is coordinatized by (u, 6, 4>,p^), where <f> is the scalar field 
value and p^ its momentum. The physical Hilbert space is simply a product 

H = Hg ® U4, , H4 = L 2 (R,d4>) , (22) 

and the Hamiltonian is Hg + H^, where is a polymer quantized scalar field 
Hamiltonian [10] . The physical evolution is given by the Schrodinger equation ^) and 
physical observables can be taken as the kinematical ones of [3], say V, <j>. Furthermore 
the mappings into auxiliary systems found for H\ can be extended to this situation, 
thus casting the physical system into one that is at least numerically manageable. 

Effective Friedman equation: Another application is a derivation of the effective 
Friedmann equation. Here we have the advantage that LQC kinematical operators 
become physical ones. Given a physical state, the derivation requires computation of 
the expectation values of the square of the Hubble operator and of the physical energy 
density operator of the dust. The latter is given by the operator po '■= — '■V~ 1 Hg ■ 
and the inverse volume operator has a suitable definition that descends from LQG. In 
a chosen factor ordering, the Hubble operator H and matter density p have the simple 
forms 

H = ni^a' 1 sin(26) , p = p c sin 2 (6) . (23) 

This allows us to express H 2 in terms of p. Taking the expectation values one arrives 
to the modified Friedmann equation 

<*>"-¥<*>> I 1 -*?) 

where an, cr p are the dispersions of H and p respectively. 

This equation is exact for all elements of T-Lg an d generalizes to any matter type. 
For if the term in (jj) is nonzero, the above equation remains valid: the formula 
for energy density is unchanged, but the density is now to be interpreted as that of 
all non-gravitational matter: p tot := — \V~ l Ha'-- The form of (|24|) implies a bounce 
for all matter for which energy density grows unboundedly with V^ 1 . 

Classical effective dynamics: Our dust frame formalism can also be applied to 
the "effective" approach developed in [TT], where it suggests a useful improvement. 
There one begins with a pair of canonical variables (x,p) (or set of pairs) subject to 
constraints, and defines quantum evolution via a set of equations of motion for the 
quantities (x), (p) and G m ' n — ((x — (x)T) m (p — (p)l) n ). Provided the constrained 
system is under sufficient control, the equations for expectation values can describe to 
a good precision the dynamics for sufficiently large classes of physical states. However 
implementation of this approach requires a procedure for converting kinematical 



87rGcrp 
3 p c 



(24) 
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observables to the physical ones. Explicit deparametrization produces a square root 
problem which requires careful approximation, and this has not yet been satisfactorily 
achieved. In the dust frame presented here, this problem vanishes allowing this 
approach to be unambiguously implemented. 

5. Summary 

We have presented applications to cosmology of a formulation of quantum gravity 
with dust time in which the gravitational part of the Hamiltonian constraint becomes 
the true physical Hamiltonian. This has several consequences which follow directly 
from the fact that the kinematical Hilbert space and observables become physical. 
Foremost among them is that the physical Hamiltonian is not a square root and so 
allows an analytical treatment in cases where it was not possible before, and permits 
a clear numerical approach where necessary. 

We studied several examples in the WDW and LQC with cosmological constant, 
all of which were shown to be analytically solvable. We also outlined possible further 
developments, including extensions to any matter type and potential, which are not 
possible with scalar field time. 
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